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Abstract 

New neutron diffraction and inelastic scattering experiments are used to 
investigate in detail the field dependence of the magnetic structure and 
low-energy spin wave spectrum of the Dzyaloshinskii-Moriya helimagnet 
Ba2CuGe207. The results suggest that the previously proposed model for 
the magnetism of this compound (an ideal sinusoidal spin spiral, stabilized 
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by isotropic exchange and Dzyaloshinskii-Moriya interactions) needs to be 

refined. Both new and previously published data can be quantitatively 

explained by taking into account the Kaplan-Shekhtman-Entin-Wohlman- 

Aharony (KSEA) term, a special magnetic anisotropy term that was predicted 

to always accompany Dzyaloshinskii-Moriya interactions in insulators. 
75.30.Et,75.10.Hk,75.30.Ds,75.30.Gw 
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I. INTRODUCTION 



The recently discovered spiral magnet Ba2CuGe207 is one of many materials known to 
have incommensurate magnetic structures.^! A fortunate combination of properties how- 
ever make Ba2CuGe20y a particularly useful model system for both experimental and 
theoretical studies of incommensurate magnetism: 1) unlike the extensively studied rare- 
earth compounds,! Ba2CuGe207 is an insulator, and thus can be conveniently described 
in terms of localized spins; 2) Helimagnetism in Ba2CuGe207 is caused by the somewhat 
exotic Dzyaloshinskii-Moriya off-diagonal exchange interactions!!! that involve only nearest- 
neighbor spins. This is in contrast with such well-known systems as Mn02 (Ref. ^) and NiBr2 
(Ref. |?]f|), where the magnetic incommensurability results from a competition between ex- 
change interactions for different neighbor pairs (geometric frustration). 3) Compared to 
such classic Dzyaloshinskii-Moriya helimagnets as MnSi (Ref. |9|,|ToD and FeGe (Ref. |lT|) , 
Ba 2 CuGe207 has a rather low (tetragonal) crystal symmetry. The result is a much richer 
field-temperature phase diagram. In particular, in Ba2CuGe207 a magnetic field applied 
along the unique tetragonal axis induces a peculiar Dzyaloshinskii-typeS incommensurate- 
to-commensurate transitionliE Applying a magnetic field in the tetragonal plane does not 
change the length of the magnetic propagation vector, but leads to its re-orientation. §4) 
The spin arrangement in Ba2CuGe207 is a perfect square lattice. This fact allowed us previ- 
ously to describe the static properties of this remarkable system using a simple and elegant 
macroscopic free energy functionalillHll 5) Last but not least, the scale of energies and wave 
numbers that characterize magnetic interactions in Ba2CuGe207 are very convenient for 
neutron scattering measurements. Magnetic fields in which the most interesting magnetic 
phase transformations occur are also readily accessible using standard equipment. 

As described in detail elsewhere (Refs. p]JT^), the principal feature of Ba 2 CuGe207 is 
a square-lattice arrangement of Cu 2+ ions in the (a, b) plane of the tetragonal non-centric 
crystal structure (space group Pk2\m, a = 8.466A, c = 5.445A). Nearest-neighbor in-plane 
antiferromagnetic exchange coupling [along the (1,1,0) direction] is by far the strongest 



3 



magnetic interaction in the system (J fa 0.96 meV per bondllS). The interaction between 
Cu-spins from adjacent planes is much weaker and ferromagnetic (J± ~ —0.026 meV 
per bond). The magnetic structure can be described as an almost- ant if erromagnetic spiral 
(Fig. 1, insert), with spins confined in the (1, 1, 0) plane and the magnetic propagation vector 
(1 + C> C> 0)j C ~ 0.0273, (1, 0, 0) being the Neel point. It was previously demonstrated that 
the helimagnetic state is stabilized by nearest-neighbor Dzyaloshinskii-Moriya interactions 
that for two interacting spins Si and S2 can be written as (S± x S2) • D^ 1 ' 2 ' . For the Cu- 
Cu bond along the (1,1,0) direction (x-axis) the Dzyaloshinskii vector D is pointing along 
(1,1,0) (y-axis), inducing a relative rotation of the interacting spins in the (x, z) plane 
[the z axis is chosen along the c -axis of the crystal]. The rotation angle (f) (relative to a 
perfect antiparallel alignment) is related to the magnetic propagation vector by <p = 2tt( fa 
0.172. Obviously, two types of domains, with equivalent propagation vectors (1 + C>C>0) 
and (1 + (, 0) will always be present in a macroscopic sample. 

By now, a large amount of experimental and theoretical work has been done on 
Ba2CuGe207, mainly dealing with the phase transitions and static magnetic properties. 
Some important issues remain unresolved however. For example, it was predicted that 
applying a magnetic field along the unique axis should give rise to a distortion of the 
ideal spiral structureJillil This so-called soliton phase is characterized by the appearance 
of higher-order magnetic Bragg harmonics. To date these additional Bargg reflections have 
not been observed directly in an experiment. As far as the spin dynamics is concerned, 
only the near-zone-boundary spin wave dispersion relations were studied. For the physics of 
the incommensurate state, it is the the low-energy, small-Q spin excitations that are most 
relevant. In the present work we continue our studies of Ba2CuGe207, investigating the field- 
dependence of higher-order magnetic Bragg peaks and the low-energy spin-wave spectrum 
in both the incommensurate and commensurate states. We find that even in the absence 
of an external magnetic field the spiral structure is distorted by the presence of previously 
disregared magnetic anisotropy. Our new results are consistent with the theoretical predic- 
tions of Kaplan,0 Shekhtman, Entin-Wohlman and AharonyBEl who demonstrated that 
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a generic anisotropy term must inevitably accompany Dzyaloshinskii-Moriya interactions. 
The new understanding of the physics of Ba 2 CuGe207 enables us to refine our interpreta- 
tion of previously obtained experimental data. A brief report on some of our new results is 
published elsewhere^. 



II. EXPERIMENTAL 

Neutron diffraction and inelastic neutron scattering measurements were performed in 
two series of experiments, on the IN-14 3-axis spectrometer at the Institut Laue Langevein 
(ILL) in Grenoble, and the SPINS spectrometer at the Cold Neutron Research Facility 
at the National Institute of Standards and Technology (NIST). Single-crystal samples of 
Ba 2 CuGe20"7 rarely survive more than one cooling to low temperatures. Two different 
crystals, prepared by the floating-zone method, were used in the two experimental runs. 
Sample A, a cylindrical single crystal of dimensions 4 x 4 x 20 mm 3 was used in experiments 
at IN-14, but spontaneously disintegrated during subsequent storage. Sample B was used 
in the second experiment on SPINS and was approximately 6 x 6 x 50 mm 3 . Crystal mosaic 
was around 0.35° FWHM for sample A and 1.2° FWHM for sample B, as measured in 
the (a, b) crystallographic plane. The mosaic spread in the perpendicular direction was 
measured for sample B and found to be around 2° FWHM. The samples were mounted 
on the spectrometers with their c -axes vertical, making (h, k, 0) wave vectors accesible for 
measurements. In both experiments the magnetic field was produced by standard split-coil 
superconducting magnets. The alignment of the c axis of the crystal with the direction of 
the magnetic field, previously shown to be crucial for high-field measurement s,0 was around 
1.4° in both runs, as measured at low temperatures. The measurements were performed 
in the field range 0-2.5 T. The sample environment was a pumped- 4 He cryostat for the 
ILL experiment and a cryopump-driven 3 He cryostat at NIST. The data were collected at 
temperatures in the range 0.35-5 K. As observed previously, cooling the sample through Tjv 
in an H « 1 T magnetic field always resulted in a single-domain magnetic structure. 
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The spin wave dispersion was measured in constant-Q scans in the range of energy 
transfers 0-0.8 meV. Neutrons of 3.5 meV or 2.5 meV fixed incident energy were used in 
most cases. Alternatively, a 3.5 meV fixed final energy setup was exploited. A Be filter was 
positioned in front of the sample to eliminate higher-order beam contamination. 40' — S — 
40' — A — 40' collimations were utilized in both runs. The typical energy resolution with 
3.5 meV incident energy neutrons was 0.075 meV FWHM, as determined from measurements 
of incoherent scattering from the sample. 

III. RESULTS 
A. Higher-order Bragg reflections 

In previous studies the only magnetic elastic peaks observed in Ba2CuGe207 were those 
corresponding to an ideal sinusoidal spiral structure with propagation vectors (1±C> ^Ci 0). 
These reflections, whose intensity seems to account for almost 100% of the expected mag- 
netic diffraction intensity appear below T/v ~ 3.2 K around antiferromagnetic zone-centers 
(h,k,0), h, /c-integer, h + k-odd. In the present study, careful elastic scans along the 
(1 + x, x, 0) line in reciprocal space revealed the presence of additional extremely weak 
peaks at (1 ± 3£, ±3£, 0), as shown in Fig. [II. These peaks are clearly of magnetic origin, as 
the temperature dependence of their intensity is similar to that of the principal magnetic 
reflections at (1±C? ±C? 0)- The additional 3rd-order peak was observed in all magnetic fields 
in the range 0-1.9 T, and always appears at (1 ± 3(, ±3£, 0), where (, defined by the posi- 
tion of the principal magnetic Bragg peak, is itself field-dependent For < H < 1.7 T 
the measured field-dependence of ( is in total agreement with previous studies. For the 
purpose of convenience we shall define Q n _ w = (1,0,0) (antiferromagnetic zone-center), and 
Qo = (C)C?0)- 111 this notation the 1st and 3rd order magnetic reflections correspond to 
momentum transfers Q n>n ± qo and Q n<n ± 3qo, respectively. 

As observed previously, at H = Hi « 1.7 T (at T — 0.35 K) the system goes through a 
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magnetic transition to a new phase that is characterized by the appearance of a new peak at 
the commensurate (1,0,0) reciprocal-space position. As discussed previously, this new phase 
may or may not be a result of the slight misalignment of the magnetic field relative to the 
c-axis of the crystal. In the present work we did not investigate this "intermediate" phase in 
detail, performing most measurements in the field ranges < H < Hi and H > H c k, 2.2 T, 
where H c is the magnetic field at which the structure becomes commensurate. 00 

In scans along the (1, 1,0) direction, shown in Fig. [I], the widths of both 1st and 3rd- 
order peaks are resolution-limited. This is not the case for transverse scans along (1, 1,0), 
where the 1st, and especially the 3rd harmonic are visibly broader than the experimental 
resolution (Fig. |^). The observed peak width pattern is consistent with both peaks having a 
zero longitudinal and a 20° transverse intrinsic widths as seen from the Q n , n reciprocal- space 
point. The transverse intrinsic Q-width of the 3rd harmonic is thus three times as large as 
that of the first harmonic. This result does not appear to depend on the applied magnetic 
field or the history of the sample. The large observed transverse width is likely to be related 
to the previously established fact that the spiral structures propagating in any direction 
in the (001) plane have almost identical energies. The (1, 1, 0) direction is only slightly 
energetically preferable.!^ Due to pinning or even thermal fluctuations, in a macroscopic 
sample an entire ensemble of spiral structures with propagation directions fanning out around 
(1, 1, 0) will therefore be realized, producing substantial transverse peak widths. 

The field dependence of the integrated peak intensities was measured in both field-cooling 
and zero-field cooling experiments. Consistent results were obtained in both types of mea- 
surement, and no signs of hysteresis were observed. In the ILL experiment the propagation 
direction of the spiral, always along (1, 1, 0) at H = 0, was found to deviate by as much 
as several degrees from this direction in higher fields. This effect is clearly due to a slight 
misalignment of the magnetic field relative to the c-axis, and the possibility to almost freely 
rotate the magnetic propagation vector in the (a, b) crystallographic plane. In the NIST 
experiment such a deviation was not observed, thanks to a slightly different and more "for- 
tunate" setting of the sample. The field dependence of the peak intensities was therefore 
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measured in this second experimental run, but, just in case, at each field, both the 1st and 
3rd-order peaks were centered in a series of transverse and longitudinal scans. The measured 
integrated intensity of the 1st and third-order reflection, as well as that of the commensurate 
peak at (1,0,0), are plotted against magnetic field applied along the c-axis in Fig. |. The 
total intensity of all three features is field-independent within experimental error. 

As seen in Fig. |, the intensity of the (1 + C, C, 0) magnetic reflection is almost field- 
independent in the range < H < H\. This appears to be in contradiction with previous 
measurements (Ref. [14], Fig. 3d), where a gradual decrease of the intensity of the 1st harmonic 
was observed with increasing magnetic field. However, we now know what was wrong with 
these previous measurements: the possibility of the propagation vector deviating from the 
(1,1,0) direction was not taken into account. In a slightly misaligned sample the field- 
induced drift of the magnetic reflections away from the line of the elastic scan was incorrectly 
interpreted as a decrease of peak intensity. Note that in the present study the centering of 
the peaks at each field ensures that this problem, even if present, does not influence the 
measurements. 

B. Spin waves 

All inelastic measurements were done in the vicinity of the (1,0,0) antiferromagnetic 
zone-center (Q « Q n ,ir)- The spin wave dispersion was measured along the (1 + e, e, 0) 
direction (x-axis). In most cases the sample was field-cooled to eliminate the need to deal 
with inelastic signal coming from the two magnetic domains. All scans at Q — Qn^ were 
repeated using zero-field-cooling to ensure that no hysteresis effects influence the measure- 
ments. The most important limiting factor in these inelastic studies is the presence of 
incoherent scattering and magnetic Bragg tails, centered at zero energy transfer. This unde- 
sirable contaminations is absent for energy transfers (hu) ^ 0.15 meV, where reliable data 
could be collected. Typical constant-Q scans obtained in the ILL and NIST experiments are 
shown in Fig. [|. Peak positions were determined by fitting Gaussian profiles to the data. 
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All the inelastic peaks studied were found to be resolution-limited. The focusing conditions 
are considerably more favorable at Q — (1 + e, e, 0), e > 0, where most of the measurements 
were performed. 

1. Zero field 

The dispersion relation measured in zero applied field is plotted in symbols in Fig. ||. One 
clearly sees three distinct branches of the spectrum. These we shall label by the wave vectors 
to which they extrapolate at zero energy transfer: Q^n ±<?o and Q n>n , correspondingly. An 
obvious and very interesting feature is the "repulsion" between the Q nt7r ±qo branches at their 
point of intersection Q = Q W)7r . Its magnitude is given by the splitting 25 ~ 0.12(1) meV. 
This effect again manifests itself at Q nj7T + 2q , where it is seen as a discontinuity in the 
Q ni7T ± qo branch. Simple empirical fits to the data (not shown) allow us to estimate the 
spin wave velocity Co ~ 5.21(3) meVA. This value is in reasonable agreement with the 
estimate c = Ja/y/2 « 5.75 meV A, obtained using the classical formula (3) in Ref. ^ and 
the exchange constant J ~ 0.96 meV, previously determined from measuring the spin wave 
bandwidth. The main characteristic of the Q nt7r branch is the energy gap ]7r « 0.18(1) meV 
at the antiferromagnetic zone-center Q n>% . 

2. Field dependence in the incommensurate phase (H < H\) 

In Figure |6] we show the spin wave dispersion measured in Ba2 CuGe2C>7 in a H = 1 T 
magnetic field applied along the c axis of the crystal. In this case the incommensurability 
parameter ((H = 1 T) = 0.0252(5). The Q n> „ ± qo dispersion curves are very similar to 
those measured in zero field and appear to be adequately described by the same spin wave 
velocity and splitting parameter 25 nj7r . Compared to the zero field case however, at H = 1 T 
the central Q n)7T branch is visibly flattened at its minimum. The gap A W)7r in this mode is 
equal to ~ 0.24 meV. Comparing this to A^ = 0.18 meV at H = 0, we find that, to a good 
approximation: 
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A^(tf) 2 = + (2g c S^ B H) 2 , (1) 

where g c = 2.474 is the c-axis diagonal component of the gyromagnetic ratio for Cu 2+ in 
Ba 2 CuGe20y (Ref. ^TJ), S — 1/2 is the spin of Cu 2+ ions and [i B is the Bohr magneton. At 
H = 1 T the measured dispersion curve for the Q n)7T branch has a new feature, namely a 
discontinuity at Qn^ + qo- This splitting, that we shall denote as 25 qo , is roughly 0.05 meV. 

In a magnetic field H = 1.5 T~ H c \ = 0.0232) the spectrum becomes substantially 
more complex (Fig. |7|). The two Q n>n ± qo modes remain essentially unchanged. The Q m>7: - 
gap in the central branch is Aq^^(H) pa 0.28 meV, which is consistent with Eq. [I]. The 
discontinuity at + qo is clearly visible in H = 1.5 T data: 25 qo ~ 0.11 meV. A new 
feature of the spectrum that is not visible in lower applied fields is the presence of a new 
excitation branch that at is seen at at (hu) ~ 0.45 meV. The shortage of beam time 
prevented us from following this branch to lower energy transfers (where its intensity should 
increase) at wave vectors where it would appear focusing: all measurements were done at 
Q = (1 + e, e, 0), e > 0. The limited data that we have at this stage is totally consistent 
with the new branch being a replica of the Q n>n mode, but centered at Q n<7r ±2qr 0) as shown 
by the corresponding solid lines in Fig. [7]. 

3. High field: commensurate phase. 

The dispersion relations measured at H = 2.5 T, well above H c pa 2.2 T, are shown 
in Fig. |^. As expected for the commensurate state, only two branches are present. Two 
peculiarities are to be noted here. First, the measured spin wave velocity Co = 4.83(3) meVA 
is significantly smaller than that seen at H < H c . Second, the gap in the higher-energy 
branch (pa 0.45 meV) is too large to be accounted for by the effect of magnetic field alone 
(2g c S^BH = 0.36 meV). If for this branch we can write: 

(huf = A 2 C + (2g c Sft B H) 2 + c\q\ (2) 
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for the "additional" gap in the commensurate phase we obtain A c = 0.28(1) meV (solid 
lines in Fig. |8|). 

IV. THEORY 

Most of the magnetic properties of Ba2CuGe207 reported to date appeared to be rather 
well described by a simple spin Hamiltonian that included only nearest-neighbor Heisenberg 
antiferromagnetic exchange interactions and the Dzyaloshinskii-Moriya cross-product terms. 
For reasons that will shortly become apparent we shall refer to this construct as the "DM- 
only" model for Ba 2 CuGe 2 7 . For a single Cu-plane in Ba 2 CuGe 2 7 the Hamiltonian takes 
the form: 

n,m ' ^n+l,m ^n,m ' ^ n,m+l) 

n,m 

+ D [(<S'n,m x Sn+l,m)y + (Sn,m x S n ,m+l ) a;] ] (3) 

Here the indexes n and m enumerate the Cu 2+ spins along the x and y axes, respectively, 
S n>m are the site spin operators, J is Heisenberg exchange constant and D is the norm of the 
Dzyaloshinskii vector. Microscopically, the Heisenberg term Ti^ represents the Anderson 
superexchange mechanism.il It arises from virtual non-spin-flop hopping of two electrons 
onto a non-occupied orbital, where they interact via Pauli's exclusion principle. As shown 
by Moriyajl the cross-product term 7i( DM ) originates from spin-flop hopping, which is made 
possible by spin-orbit interactions. 

The classical ground state of the DM-only model is an ideal sinusoidal spin-spiral. The 
experimental observation of higher-order magnetic reflections in zero magnetic field tells us 
that this model is not a fully adequate description of Ba 2 CuGe 2 07: something is missing from 
the Hamiltonian ([3]). To understand what is going on we first note that for two spins S\ and 
S2, interacting via isotropic exchange and the Dzyaloshinskii-Moriya term, the interaction 
energy is minimized at — \J J 2 + D 2 S 2 when both spins Si and S2 are perpendicular to D, 
forming the angle 7r + a, where a = arctan (.D/J). Therefore, the Dzyaloshinskii-Moriya 
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cross-product term 7i( DM ) lifts the local 0(3) symmetry of the Heisenberg Hamiltonian and 



creates an effective easy plane anisotropy of strength a/ J 2 + D 2 — J ~ D 2 /2J. 

Relatively recently Kaplan0 and Shekhtman, Entin-Wohlman, and Aharonylli (KSEA) 
argued that in most realizations of Moriya's superexchange mechanism this apparent easy- 
plane anisotropy is an artifact of the omission of terms quadratic in D in the expansion of 
the true Hamiltonian of the system. If such terms are properly included, the 0(3) symmetry 
of a single bond is restored by an additional term (V J 2 + D 2 — J)/ ' D 2 (Si ■ £>)(/*>2 ■ D) ~ 
(1/2 J) (/Si ■ D)(S2 ■ D)- Note that this additional interaction has the form of easy-axis 
two-ion anisotropy and its strength is such that it exactly compensates the easy-plane effect 
of the Dzyaloshinskii-Moriya cross-product. With this term included, the ground state of 
the two-spin Hamiltonian has full 0(3) symmetry. The energy of two interacting spins 
pointing parallel and antiparallel to D, respectively, is exactly equal to that of two spins 
perpendicular to D and forming the angle ir + a between themselves. We shall refer to this 
"hidden symmetry" term as the KSEA anisotropy term or KSEA interaction. For a recent 



discussion of this subject see Ref. 23 



To properly account for KSEA interactions in our model of Ba2CuGe207, where D ^ J, 
the spin Hamiltonian can be rewritten as follows: 

ft = ft(H) + ft(DM) + ^(KSEA) = 

^ ; [(/ (S n m • S n ^i m -\- S n<m • S n,m+l) D ({<Syj jJ7l X S n ^.\ m )y -\- (S n<m X S ' n,m+l) x) 

(4) 



n,m 
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D 1 

j ( cy qy I qx qx > 

c2j\ u n,m u n+l,m ' u n,m n,m-\-lj 



Can this Hamiltonian (the "DM+KSEA" model) account for both new and previously pub- 
lished experimental data on Ba2CuGe20y? In the following sections we shall systematically 
investigate the effect of the KSEA term on static and dynamic properties of a DM helimag- 
net, and show that indeed it can. 



A. Static properties 
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1. Free energy in the continuous limit 



As the period of the spiral structure in Ba 2 CuGe 2 7 is rather long (« 36 lattice spacings), 
we can safely use the continuous approximation to describe it0El In this framework the 
magnetic free energy is expanded as a functional of a slowly rotating unitary vector field 
n(r). At each point in space n(r) is chosen along the local staggered magnetization. The 
Hamiltonian (Rf) then gives rise to the following free energy functional: 



F(° M ) = / dxdy [| {{d x n -^-e y xn) 2 + (d y n - | e x x nf 
a 2 2 , (X± ~ X\\)(H ■ n) 2 X ±H 



(5) 



2 A 2 ' 

In this formula p s is the spin stiffness, that in the classical model at T = is given by 



p s = S 2 ^/ J 2 + D 2 W S 2 J, a as before is the equilibrium angle between two spins defined 
as a = arctan(D/ J), A is the nearest-neighbor Cu-Cu distance, x\\ an d X± are ^ ne local 
longitudinal and transverse magnetic susceptibilities, respectively. Their classical T = 
values are x± — (5 , / i s) 2 /(4^A 2 ) and x\\ — 0, correspondingly. In Eq. |5] we have included the 
Zeeman term that represents interaction of the system with an external magnetic field H . 

The term — a 2 p s n 2 / 2A 2 in Eq. ^ deserves some comment. It has the form of a magnetic 
easy-z-axis anisotropy and represents the combined effect of the effective (xz) and (yz) easy 
planes produced by DM interactions on the y- and x-bonds, respectively. This term is 
eliminated by KSEA interactions that modify Eq. || as follows: 



F (DM + KSEA) , dxdy 

, (x±-X\\)(H-n) 2 X ±H 2 



y (^(d x n -^e y xn) 2 + (d y n ('-.,■ x nf ) 



(6) 



2 2 

This equation is in agreement with Eq.(3) in Ref. [15|. Comparing Eqs. |^ and |] one concludes 
that in the continuous limit for the square-lattice spin arrangement found in Ba2CuGe207, 
KSEA interactions (two sets of easy axes, for x and ?/-bonds, respectively) are indistinguish- 
able from an overall easy-(x?/)-plane anisotropy of relative strength 5 = a 2 /2. 
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In this work we are mostly concerned with the effect of a magnetic field applied along the 
[001] crystallographic direction, i.e., along the z-axis. Under these conditions the propaga- 
tion direction of the spin spiral in Ba2CuGe207 is either along the x or y axis (two domain 
types are possible). Moreover, as we shall prove rigorously while discussing the spin waves 
in the system, the magnetic structure remains planar despite the two types of Dzyaloshin- 
skii vectors, along the x and y axes (for the y and x bonds, respectively). This fact allows 
us to write the components of vector n(r) as (sin#(x), 0, cos#(x)), where 6(x) is the angle 
between local staggered moment n(r) and the z axis, for a helix propagating along the x 
direction. The free energy can be then rewritten in terms of the 9{x) as 



F (DM+KSEA) = J dxdy 



Ps {d x 9-{a/K)f ,a 2 p s ( x± - Xll )H\ 2 X ±H" 
2 + ( 2A^ + 2 )C ° S — 



• (7) 



This is exactly Eq. (1) of Ref. ^ modified to include the effects of an easy (xy) plane 
anisotropy p s a 2 n 2 z /2h 2 = const — p s a 2 cos 2 9 /2A 2 , coming from the KSEA interaction on y 
bonds. As seen from this equation the sole effect of such anisotropy is to renormalize the 
external field to 



H eS (H) = J IP + a 2 p s /A 2 ( X ± - (8) 



2. Critical field and magnetic propagation vector 

One important consequence of what is said above is that all our previous results, obtained 
in Refs. |i~3| , |i~4T , can be recycled by substituting H e s(H) for H in all formulas. Our conclu- 



sions regarding the field-induced commensurate-incommensurate Dzyaloshinskii transition 
in Ba2CuGe207 remain valid in the presence of KSEA interaction. The KSEA interaction, 
however, modifies the value of the critical field H c . Indeed, substituting H ea {H) for H in 



the Eq. 5 of Ref. one gets ^# c 2 + a 2 /A 2 p s /( X ± - X\\) = (™/^)\J Ps/(x± ~X\\)- /From 
this we immediately obtain: 



H c = a^B x p^. (9) 

2A V X± ~ XII 1 ^ 
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We see that the KSEA term reduces the critical field by the universal factor Jl — A/n 2 ~ 
0.771. 

In order to obtain the field dependence of the inverse period of the structure ( one has 
to rewrite Eqs. (4,7) of Ref. |l3| as 

2ir((H) ir 2 



a 4E{/3)K(J3) 
Hes(H) f3 



(10) 
(11) 



H cS (H c ) E{{3) 

Here f3 is an implicit variable. In case when the deformation of the spiral is weak ((a — 
2n((H))/a <C 1) one can safely use the linearized formula: 

— — — - = 1 — — ( 7777—777-7 ) + higher order terms. (12) 
a 32 \2H eS (H c ) J 

From this formula one can derive by how much KSEA interactions deform the spiral in zero 

external field. Let us define <fi = 2n((0) as the average angle between spins in the spiral in 

zero field. This parameter is easily accessible experimentally and equal to <fi = 2tt ■ 0.0273 = 

0.172 ~ 10°. 

Recalling that H ef[ (0)/H eS (H c ) = 2/n and plugging it in Eq. [12] one gets <p/a ~ 1 — 1/32 

or 

32 32 

a = arctanp/J) = -0 = — 2ttC(0). (13) 
The KSEA term thus increases the period of the structure in zero field by roughly 3%. 

3. Higher- order Bragg harmonics 

An important implication of Eq. § is that even in zero applied field the effective field is 
non-zero. The result is that the spiral structure is distorted even in zero field and higher- 
order (odd) magnetic Bragg peaks are present. To obtain a theoretical form for the field 



dependence of the 3rd harmonic we can use Eqs. 17, 18 in Ref. [T3|. More practical than 
the resulting expression is its linearized form, that applies in the limit (a — 2tt((H)) a 
(weakly distorted spiral) : 
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1 ( kH cS {H) 
h 256 \2H cS (H c ) 



256 



(14) 



Here Ii and /3 are the intensities of the first and third harmonic, respectively. One can see 
that for H = the third harmonic is predicted to be smaller than the first one by a factor 
of 1/256 ~ 3.9 x 1(T 3 . 

Comparing Eqs. |TJ] and one can see that for weak distortions the intensity of the third 
harmonic is proportional to the relative decrease of ({H): 



CJH) 

C(o) 



l - 



8(I 3 (H) - I 3 (0)) 
h 



(15) 



4- Comparison with experiment 

We now have to make sure that our results for the DM+KSEA model are consistent 
with both the previously published (Ref. 0,0) and new neutron diffraction results on 
Ba 2 CuGe 2 7 . First, let us compare the previously measured ((H) curve with the predictions 
of Eqs. 0,0. The experimental data for the incommensurability parameter (Ref. [13]) is 
plotted against if 2 in Fig. []. The solid line is the prediction of Eqs. |TU| , |rT| . The dashed 
line is the theoretical curve previously obtained without including the KSEA term in the 
Hamiltonian.@0 In plotting both these curves we have assumed the actual (measured) 
values for H c = 2.15 T and £(0) = 0.0273. Within experimental statistics it is practically 
impossible to distinguish between the two theoretical dependences and the data fits both of 
them reasonably well. 

While it appears that the shape of the ((H) curve can not be used to extract information 
on KSEA interactions, the actual numerical value of H c in the DM-only and DM+KSEA 
models is substantially different. For the low-temperature limit in Ba2CuGe20y we can use 
the classical expressions p s = JS 2 = 0.24 meV, x\\ — and x± = (fi'cA i B) 2 /8JA, where 
g c = 2.47 is the c-axis gyromagnetic ratio for Cu 2+ in Ba2CuGe2O7.0 One can expect these 
classical estimates to be rather accurate, as they rely on the effective exchange constant 
J, that itself was determined from fitting the classical spin wave dispersion relations to 
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inelastic neutron scattering data.El Substituting these values into Eq. ||we immediately obtain 
^dm+ksea) — 2 T. This is much closer to the experimental value H c = 2.15 T, than the 
estimate Hp M ^ = 2.6 T for the DM-only model. 

The new data for the field dependence of relative intensities of the first and third Bragg 
harmonics becomes consistent with theory only if KSEA interactions are properly taken 
into account. Indeed, the KSEA term is necessary to reproduce the observed distortion of 
the spiral in zero magnetic field. In Fig. 0(b) the solid lines are plotted using Eq.|l4| and 
H c = 2.15 T. The dashed lines are results for the intensity of the third Bragg harmonic ob- 
tained previously for the DM-only model.0 Clearly the DM+KSEA model gives an excellent 
agreement with experiment (open symbols in Fig. 0(b)), while the DM-only Hamiltonian 
fails entirely to account for the available data. 

In Fig. 0(c) we check the validity of theoretical prediction of Eq.|15|, which is supposed to 
hold both with and without KSEA terms. The excellent agreement of theory and experiment 
confirms the validity of our picture of weakly deformed almost sinusoidal spiral. 



B. Spin dynamics 

We now turn to calculating the classical spin wave spectrum in the DM+KSEA model 
for Ba 2 CuGe207. This task will be accomplished in several separate steps. First, we shall 
derive the spectrum for a square-lattice Heisenberg Hamiltonian, including DM interactions 
only for the x-axis bonds. Second, we shall consider the effect of DM interactions along 
the y -axis bonds, showing that they do not disturb the planar spiral structure and do not 
influence the dispersion relation along the x-direction. Next we shall analyze the effect of 



adding the KSEA term, following the method described in Ref. ^J,^. While at this stage 
we do not have results for spin wave dispersion in the DM+KSEA model in the presence 
of an arbitrary external magnetic field, we shall consider the case H > H c and derive an 
expression for A c - the anisotropy gap in the commensurate state. 
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1. Dzyaloshinskii-Moriya interactions for x-axis bonds only. 



We start from a truncated version of the Hamiltonian (|j): 

)vl (16) 



n,m 



It is easy to see that classically this Hamiltonian is minimized by a perfect helicoid propa- 
gating along the x-axis with all spins lying in the (xz) plane: 

(S z ntm ) = (-l) n+m (S)cosna; 

(Sl m ) = (-l) n+m (S)smna. (17) 

The standard procedure to calculate the spin-wave spectrum is to rewrite this Hamilto- 
nian in terms of spin projections on the new rotating coordinate system, where the direc- 
tion of the equilibrium value of spin at (n, m) defines the local z' axis in such a way that 
(Sum) = (— l) n+m (5')- We leave the y coordinate unchanged, and select the new x' axis 
to be orthogonal to both z' and y' = y. Substituting S* m = S* m cosna — S^' m smna, 

Sn,m = Sn,m cosna + Snm s * nna > an< l ^n,m = ^n,m m the Hamiltonian flip] ) and using 
a = arctan(.D/J) we obtain: 

^ ■* = X] W^ 2 + D 2 (S nm S n+l m + S n ^ m S n+1 ^ m ) + JS^ m Sn+i ;m + JS n m ■ S n ^ m+1 ]. (18) 

n,m 

In these coordinates the Hamiltonian is simply that of a square lattice AFM with easy-plane 
exchange anisotropy on bonds along the x direction. In agreement with the discussion in 
Section [IV A 1| the relative strength of this anisotropy is given by: 

>,2 



V J 2 + D 2 -J D 2 or 



(19) 



J 2J 2 2 ' 

Applying the Holstein-Primakoff formalism we write the spin projection operators as 



S( m = (-l) n+m (S-< m a„, m ), = (-ir+ m JS/2(a n , m +alJ, S« = iJS/2 (4 >r 
From Eq. flT6| ) it is then straightforward to extract the quadratic part of the spin wave 
Hamiltonian: 
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7"^ — JS (^ a n jr n a n,rn a n,m a n+l,m a n,m a n,m+l a n,m a n+l,m a n,m a n,m+l) 
n,m 

+ a n,m) { a li+l,m 

(20) 

After performing Fourier and Bogolyubov transformations to diagonalize this Hamiltonian, 
one readily obtains the spin wave spectrum: 



S(k x ,k y ) = JS^J[4: + 25 -5 cos k x ] 2 - [(2 + 5) cosk x + 2cosk y } 2 . (21) 

This spectrum has one Goldstone branch at k x = k y = 0, that corresponds to the continuous 
symmetry of a simultaneous rotation of all spins in the (xz) plane. At k x = k y = ir the 
spectrum has a finite gap AJS\f5 = 2\J~2DS due to the easy (xz) plane anisotropy coming 
from the DM without KSEA correction. 

Now we have to recall that in the above derivation the wave vectors k x , k y correspond to 
a rotating system of coordinates. They are thus distinct from the actual component of the 
scattering vector in a neutron experiment. To get the proper spin wave spectrum one has 
to perform a reverse coordinate transformation to the laboratory system: 

*~Vi m ~ m COS TICX S n m sin TlOt = 



= (-l) n+m [(S - a^ m a„ im ) cosna - ^fs/2 (a n , m + a^ m ) sinraa]; 
S* m = S*[ m cosna + s£ m sin na = 



(-l) n+m [(S - a ] n m a ntm ) sinna + \fs~/2 (a n , m + a^ m ) cosna]; 



^n,m ~ ^n,m ~ £y$/2 (ajj jTn 0>n,m)- (22) 

The x-axis dispersion of three spin wave branches in laboratory system is shown in |T0"|(a). 
The dynamic structure factor S yy (Q,uj) has a single magnon peak at the energy given by 
£(Qx->Qy) (the Q^-branch). The structure factors S xx (Q,u) and S zz (Q,u) each contain 
two magnon branches with dispersion relations given by £(Q X + vr + a, Q y + ir), and S(Q X + 
ir — a,Q y + 77 ) (the Q n>n ± go-branches). As expected, the zeroes of energy in these two 
modes are precisely at the positions of magnetic Bragg peaks at Q n>n ±q . A curious feature 
of this plot is that all three branches are nearly degenerate at the AFM zone center. 

19 



2. Dzyaloshinskii-Moriya interactions along y-axis bonds 

Let us now consider Dzyaloshinskii-Moriya interactions for the bonds in the ^-direction. 
Their contribution to the spin Hamiltonian can be written as: 

o/(2) _ ST< r)( qy qz _ q z qv \_Y^nra (qz — Q z \ - 

' L / j ^ \ kJ n,m kJ n,m+l LJ n,m LJ n,m+l) / j ±y "n,.ml IJ n.m+1 u n,m— l! 

n,m n,m 

= iDS'^2(—l) n+m sinna [a„ im a nim+ i — alim a n,m+i\ + third order terms (23) 

n,m 

The absence of terms of the first order in a„ jm and o) n m means that in the original (flat-spiral) 
spin configuration the force acting on each spin produced by the added Dzyaloshinskii-Moriya 
coupling on the y bonds is equal to zero. Thus, switching on the ?/-axis DM interactions does 
not disturb the planar helimagnetic ground state of the Hamiltonian TiP~\ which therefore is 
also the ground state of H.^ = Ti.^ + 7i( DM ) . This a posteriori verifies our assumption that 
spins continue to lie in the x — z plane in the presence of Dzyaloshinskii-Moriya interactions 
on y-bonds, made in the section (IV A 1). 

New terms quadratic in a n ^ m and a\ im are indeed introduced by the Dzyaloshinskii- 
Moriya interactions on y-bonds, and the spin wave spectrum is thus altered. After Fourier 
transformation Eq. ( p3[ ) becomes: 

iDS 

= ^ sin k y [o){k x , k y )a\—k x + tt + a, —k y + tt) — 

— a\k X) k y )a){—k x + tt — a, —k y + tt) + 
+ a(k x , k y )a(—k x + tt — a, —k y + tt) — 

— a(k x , ky)a(—k x + tt + a, —k y + tt)]. (24) 

The analysis of this term for general k y is rather complicated and should be done by matrix 
diagonalization similar to that described in the next subsection for calculating the effects 
of KSEA interactions. Fortunately, for spin waves propagating along the x-axis (k x = or 
k x = tt) the contribution of HS 2 ^ is exactly zero, thanks to the sin k y prefactor in Eq. |23j. 
In other words, as long as we are concerned with spin waves propagating along the (110) 
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direction in Ba2CuGe207 we can totally disregard the contribution of Dzyaloshinskii-Moriya 
interactions along the y-axis bonds. 



3. Influence of KSEA interactions. 

Having understood the spectrum for the DM-only model, we can proceed to include 
KSEA terms in our calculations. We first note that if our system were 1- dimensional, the 
inclusion of the KSEA term would fully restore 0(3) symmetry, making the commensurate 
and spiral phases degenerate. In terms of spin waves this would signify a complete softening 
of the magnon branch at the AF zone-center Q^. As will be demonstrated below, in 
the case of a 2-dimensional spin arrangement in Ba2CuGe20"7 the magnon softening at Q n)7T 
produced by KSEA interactions is incomplete. 

As discussed in Section [IV A3| , in the presence of the KSEA term the ground state is a 



distorted flat spin-spiral. In this situation the transition to a uniformly rotating coordinate 
system used in Section [IV B 1| loses its usefulness. Instead, we must rotate the coordinate 



system for spin quantization at each site in such a way, that the z-axis follows the rotation 
of the spins in the distorted helix: 



s n, m = (-l) n+m [(S - 4 im a n>m ) cos# n>m - (a n , m + a+J sin0„ iTO ]; 



s n,m = (-l) n+m [(S - 4 )m a„ im ) sin 6^™ + (a n>m + o) n m ) cos0 n)OT ]; 

Sn,m = Sn,m = * \f^/^( a n,m ~ a n,m)- (25) 

Here 8 n m denotes the angle between the local spin axis and z -axis in the x — z plane. The 
Hamiltonian 7v- x > + (as explained above, HS 2 ^ is not relevant to the dispersion along 

the x-axis that we are interested in) is then rewritten as: 



cos(6 n+ltm - 9 n>m - a) 



cos a 



2 ( a n,m a n.m) ( a n+l,m a ra+l,m)^ ^ ( a n,m a n,m) ( a n+l,m a n+l,in) 
— COs(# njm +i — 6 nm ) (j^S — 0^ nm a nm )(S — Q'n j . m +lO'n,m+l) + 
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~l~ 2 ^- a ri,m "I" a n,m) ( a n,m+l "I" a n,m+l)^ 

— ~^{ a n,m ~ a n,m) ( a n,m+l — a n,m+l) 

' ■ / S . j \/t "\ 

— I ^"v "",™ — a n,m) ( c *'n+l,m — a n+l,mJ + 

,m+i) sin $ n>m sin 9 n ^ m +i + 
+ (S - an, m o n , m )(5' - al m+1 a njm+1 ) cos0 n , m cos# nim+ i)| + linear terms (26) 



If all angles B n ,m in the above expression are given by the solution of the sin-Gordon 
equation determining the ground state, the linear terms will vanish: they represent a static 
uncompensated force acting on the spins and must not be present in an equilibrium spin 
configuration. In general, Eq. ^6] can not be diagonalized analytically Fortunately, we 
are dealing with a rather weakly distorted structure and can safely restrict ourselves to 
calculating the effect of the KSEA term to the first order in 5. It is easy to show that 
the easy x — y plane anisotropy of strength 5 deforms the spiral in such a way that is 
9n,m = qn + (5 / '4a 2 ) sin 2qn+0(5 2 ), where q = a — 0(5 2 ). Therefore, within our accuracy one 
can assume 9 n<m ~ an in Eq. (|26|) . In particular the anisotropy dependence of cos(6 , n+l m — 
9n,m — ot) ~ cos(g — a + (5/AJa 2 ) sm2qn + 0(5 2 )) ~ 1 + 0(5 2 ) can be disregarded. With 
these simplifications and after Fourier transformation Eq. (E26|) becomes 



4 + 5- -cos/c y j a(k x ,ky) ] a(k x ,k y )- 



( S \ a(-k x ,-ky)a(k x ,k y ) + a(-k x ,-k y Ya(k x ,k y y 
— 12 cos k x + 2 cos k y + - cos k y I — h 

+ (A + ^ cosfc^ (a(k x + 2a, k y ) ] a(k x , k y ) + a(k x - 2a, k y ) ] a(k x , ky)} + 

5 I a{—k x + 2a, k y )a{k xi k y ) + a(—k x + 2a, k y Ya(k x , k y Y 

+ — cos ky 



+ 



4 y V 2 
a(-k x - 2a, k y )a(k x , k y ) + a(-k x - 2a, k y Ya(k x , k y ) r 



2 

From this equation we can already qualitatively understand the role of KSEA interactions. 
Their main impact is the introduction of terms that couple magnons with wave vectors that 
differ by 2q . This coupling will have the largest effect when acting on a pair of magnons 
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(27) 



of equal energies. The result will be discontinuities in the magnon branches at certain wave 
vectors, that for the distorted helix become new zone-boundaries. This picture is very similar 
to the formation of a zone structure and zone-boundary energy gaps in a free electron gas, 
subject to a weak periodic external potential. Another consequence of KSEA interaction is 
the reduction of the energy gap in the Q n>ir branch at k x = k y = 7r from 2\/2DS to 2DS. 
However, contrary to the one- dimensional case this gap does not become zero, i.e. a new 
Goldstone excitation does not appear in two dimensions. This can be derived by looking at 
the part of the Eq. (|27|) , which does not involve mixing of branches separated by 2q , and, 
therefore, yields to the standard analytical calculation. 

To actually calculate the spin wave spectrum we have to find a transformation of Bose 



operators that would diagonalize the Hamiltonian (p7|). This transformation must respect 
Bose commutation relations, and for a rather general case of helimagnetic structures is 
described in detail in Refs. ^5|. It is essentially a Bogolyubov transformation involving 
a column vector of four operators: a(k x ,k y ) = (a(k x — a, k y Y , a(— k x + a, —k y ),a(k x + 
a, k y y, a{— k x — a,—k y )). The relevant part of the Hamiltonian (p7|) can be written as 
Ji = (1/2) J2kx,k y ^{kx, k y YVa(k x , k y ) , where the 4x4 matrix V is given by 



/ 



V 



\ 



\ 
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A(k x - a, k y ) B(k x -a,k y ) C(k x -a,k y ) D(k x -a,k y ) 

B(k x -a,k y ) A(-k x + a,-k y ) D(k x -a,k y ) C(-k x + a,-k y 

C(k x + a,k y ) D(k x + a,k y ) A(k x + a,k y ) B(k x + a,k y ) 

D(k x + a,k y ) C(-k x - a,-k y ) B(k x + a,k y ) A(-k x - a,-k y 

Here A(k x ,k y ) = JS [2 + 35 - (8/2)cosk y ], B(k x ,k y ) = -JS [(2 + 25)008^ + 2cosfc s/ + 
8/2 cos k y ] , C(k x ,k y ) = JS [8/2 + (5/4) cos A;^], and D(k x ,k y ) = JS (8/4) cos k y . To diago- 
nalize the spin wave Hamiltonian and at the same time ensure the conservation of commu- 
tation relations we have to find a matrix Q such that Q^VQ is diagonal, while Q^gQ = g, 
where g is the diagonal matrix with diagonal elements (1, —1,1,-1). This is equivalent to 
diagonalizing the matrix (^.0 

To obtain numerical results that could be directly compared to our measurements on 
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Ba 2 CuGe20y we used the independently measured values for J = 0.96 meV and D/J m 
arctan(D/J) = a — ||</> = 0.177. A numerical diagonalization of gV(0,0) was performed 
using Mathematica software package to yield the eigenvalues £ = 0.172 meV, £\ = 0.297 
meV and £ 2 = 0.171 meV. These are the energies of the Q n>1T (£o), an d Qn,* ± Qo 
£2) branches at the AFM zone center Q n)7T . The splitting was predicted to be 25 = 
(£1 — £ 2 ) = 0.12 meV. This value is indistinguishable from the actual splitting observed 
in Ba2CuGe2C>7, quoted in the previous section. We can also calculate the splitting in the 
Q T)7r branch at Q nt7r ± q: 25 qo = 0.049 meV. Experimentally, this splitting was not observed 
in zero field, but is small enough to be well within the experimental error bars. At higher 
fields the discontinuity at this wave vector becomes apparent (see Section ( [IIIB 2|) ). The 



gap in the Q T)7r -branch, A n ^ = £ = 0.172 meV, is also in very good agreement with the 
INS measurements. Entire dispersion branches calculated numerically using the technique 
described above are shown in Fig. 10(b). They can be also seen as solid lines in Fig. [| and 
apparently are in very good agreement with experimental data. 

4- Spin-wave spectrum in the spin- flop phase 

As we have already mentioned, we presently do not have theoretical results for the spin- 
wave dispersion in the presence of an external magnetic field. However, we can make some 
predictions for the spin- wave spectrum in the spin-flop phase (i.e., for H > H c ). After 
some tedious calculations that are omitted here, but are very similar to those performed 



in Section |IV B 2\ , one arrives at the result that the contribution of the Dzyaloshinskii- 
Moriya term for the x (y) bonds is proportional to sin k x (sin k y ) and therefore exactly 
vanishes at the AFM zone-center. In order to calculate the additional energy gap A c in 
the spin-flop phase we thus need to consider only the KSEA terms. At Q = Q T)7r (long- 
wavelength limit) the effect of KSEA interactions is identical to that of conventional easy- 
plane exchange anisotropy. The spectrum of a Heisenberg AFM with such anisotropy in 
a magnetic field is well-known.il Both field and anisotropy split the two-fold degenerate 



24 



magnons in a Heisenberg system to give a gapless mode with linear dispersion and an 
"optical" mode with the energy gap at Q n>1T given by: 



A c = 2V2JSV25 = 2V2DS. 



Substituting the known numerical values into this formula we obtain A c = 0.24 meV, which 
is in reasonable agreement with the experimental result A c = 0.28 meV. 

V. DISCUSSION 

We see that both the static and dynamic properties of Ba2CuGe20"7 are quantitatively 
consistent with the presence of KSEA interactions. To be more precise, the experimental 
data unambiguously indicate the presence of an easy-plane anisotropy of exactly the same 
strength as predicted by the KSEA mechanism. It is important to stress that in a slowly 
rotating helix it is impossible to distinguish experimentally between single-ion easy-plane 
anisotropy of type En.m {S z n r^f /2, two- ion anisotropy (£ n ,m [^, m ^+i, m + S* m 5£ iTO+1 ] /2 
) or KSEA-type anisotropy ( J2n,m &n m^n+i m + &n m^n m+i ) of the same strength. Indeed, 
the difference between a pair of easy axes (KSEA term) and an easy plane (conventional 
single-ion or two-ion anisotropy) becomes apparent only when the period of the structure is 
comparable to the nearest-neighbor spin-spin separation, i.e., is only manifested in lattice 
effects. Alternatively, it can be observed in strong magnetic fields when the canting of 
spins towards the field direction becomes substantial. In Ba2CuGe20"7 , where the magnetic 
structure has a rather long periodicity, and where even at the critical field the uniform 
magnetization (spin canting) is small, these effects are expected to be insignificant. It is 
entirely possible that the weak "quadrupolar" in-plane anisotropy seen in horizontal-field 
experiments^ is in fact such a lattice effect. Note that its strength is extremely small, of the 
order of 7 x 10 -9 eV, and yet it can be reliably measured in a diffraction experiment where a 
magnetic field is applied in the (ab) crystallographic plane at different angles to the a axisS. 
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Another possible manifestation of lattice effects is the intermediate phase seen just before 
the CI transition in a magnetic field applied along the c crystallographic axis0. A further 
study of these phenomena that distinguish between KSEA and other type of anisotropy in 
Ba2CuGe20y is an interesting topic for future experimental and theoretical work. 

One final comment has to be made in reference to another possible player in the spin 
Hamiltonian for Ba2CuGe207: dipolar interactions. In principle, the ever-present dipolar 
term can and will influence both the ground state spin configuration and the spin wave 
spectrum in Ba2CuGe207. Its effect however is expected to be insignificant compared even 
to the weakest terms that we have considered in our treatment. Indeed, nearest-neighbor 
spins in Ba 2 CuGe 2 07 are separated by A = a/y/2 ~ 6 A. For two nearest-neighbor spins 
the energy of dipolar coupling is of the order of (g/ig) 2 /A 3 pa 1 /xeV. This is an order of 
magnitude less than the smallest energy scale in our model, which is the strength of KSEA 
anisotropy Ja 2 /2 pa 15/xeV. Moreover, in an almost-antiferromagnetic structure long-range 
dipolar interactions will be heavily suppressed by the sign-alternation in the contribution of 
individual pairs of interacting spins. Our neglecting dipolar interactions in all the derivations 
above is thus justified. 

In summary we have demonstrated that KSEA interactions can result in very interesting 
measurable effects, and that no additional anisotropy is needed to reproduce the behavior 
observed in Ba 2 CuGe 2 07. 
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FIGURES 

FIG. 1. Typical elastic scans along the (1,1,0) direction in the vicinity of the antiferromagnetic 
zone-center (1, 0, 0), measured in Ba2CuGe20"7 at T = 0.35 K in zero field (top) and in a H = 1.6 T 
magnetic field applied along the (0,0,1) direction (bottom). Note the logarithmic scale on the 
y-axis. The solid lines are guides for the eye. The arrows show the positions of the principal 
magnetic Bragg reflections at (1 + £,£,0), characteristic of a spiral spin structure, and the 3rd 
harmonic at (1 + 3£, 3£, 0), a signature of a slight distortion of the helicoid. Insert: a schematic of 
the magentic structure showing a single Cu-plane in Ba2CuGe20y . 

FIG. 2. Transverse elastic scans through the lst-order (a) and 3rd-order(b) magnetic Bragg 
reflections measured in Ba2CuGe207 at T = 1.3 K and H = 1.6 T. Solid lines are Gaussian fits to 
the data. The shaded Gaussian represents the calculated experimental Q-resolution. The intrinsic 
angular width of both peaks is « 20° as seen from the (1, 0, 0) antiferromagnetic zone-center. 

FIG. 3. (b) Measured field dependence of the magnetic Bragg peak intensities in Ba2CuGe207 
at T = 0.35 K. Solid and open circles show the beha vior of the lst-order and 3rd-order incom- 
mensurate Bragg peak integrated intensities, respectively. The intensity of the commensurate peak 
at the antiferromagnetic zone-center is plotted in open triangles, (b) Square root of the ratio of 
intensities of the 3rd and 1st harmonic plotted against the square of the applied magnetic field. 
The solid and dashed lines show the theoretical prediction for the DM-only and DM+KSEA mod- 
els (c) Measured intensity of the 3rd harmonic plotted against the normalized incommensurability 
parameter (. The solid line shows the prediction of the DM+KSEA model. 

FIG. 4. Typical inelastic scans measured in Ba2CuGe207 in the two experimental runs, at 
ILL (top) and NIST (bottom), respectively. The heavy solid line is a multiple-Gaussian fit to the 
data, and the shaded curves represent the individual Gaussians. The gray area in the top panel 
shows the position of a "Bragg-tail" spurious peak. 
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FIG. 5. Spin wave dispersion curves measured in Ba2CuGe20"7 in zero magnetic field. The 
data collected at T = 0.35 K and T = 1.5 K are combined in this plot. The solid lines are 
parameter- free theoretical curves as described in the text. Dashed lines are guides for the eye and 
the solid circles on the abscissa show the positions of the observed magnetic Bragg peaks. 

FIG. 6. Spin wave dispersion curves measured in Ba2CuGe207 in a H = 1 T magnetic field 
applied along the (001) direction at T = 0.35 K. Solid lines are guides for the eye. Dashed lines 
are as in Fig. [|. 

FIG. 7. Spin wave dispersion curves measured in Ba2CuGe207 in a H = 1.5 T magnetic field 
applied along the (001) direction at T = 0.35 K. The lines and symbols as in previous figures. Note 
the additional branch in the spectrum. 

FIG. 8. Spin wave dispersion curves measured in Ba2CuGe207 in a H = 2.5 T magnetic field 
applied along the (001) direction at T = 0.35 K (commensurate spin- flop phase). The upper solid 
line is a fit to Eq. ^. The lower line is a linear fit. 

FIG. 9. Field dependence of the incommensurability parameter as previously measured in 
Ba2CuGe207. The solid curve is plotted using Eqs. [l0|JTl| , that takes into account KSEA anisotropy. 
The dashed curve is the prediction of the DM-only model. 

FIG. 10. Theoretical predictions for the spin wave dispersion along the x-axis in (a) the 
DM-only model, (b) DM+KSEA model. The effect of the KSEA term is to couple magnons 
separated by 2 go, which leads to the appearance of new gaps in the spectrum, and reduce the gap 
in branch. 
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